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INVESTIGATION OF THE MARTIAN ENVIRONMENT 
BY INFRARED SPECTROSCOPY ON MARINER 9 


ABSTRACT 

The infrared spectroscopy experiment on Mariner 9 pro- 
vides extensive information on the Martian environment, 
including spatial, diurnal, and secular dependences of 
atmospheric and surface parameters. Measurements ob- 
tained during and after the planet-wide dust storm indicate 
that large diurnal variations in atmospheric temperature 
existed up to at least 30 km; winds inferred from the tem- 
perature fields show a strong tidal component and signifi- 
cant ageostrophic behavior. With the dissipation of the 
dust, the daily maximum in the atmospheric temperature 
field moved from approximately latitude -60°and late 
afternoon local time to near the subsolar point in latitude 
and time. Analysis of spectral features due to the atmos- 
pheric dust indicates aSiC >2 content of 60 ±10%, implying 
that substantial geochemical differentiation has occurred. 
Water vapor estimates indicate abundances of 10-20 pre- 
cipitable micrometers, less than has been inferred by 
ground-based methods in similar phases of previous 
Martian seasons. Between November 1971 and April 1972 
no gross latitudinal or temporal dependence in the water 
vapor distribution has been detected from the south polar 
region to the equator. Water vapor has not been detected 
over the north polar regions. Surface pressure mapping 
has been carried out from which topographic relief of 
nearly two pressure scale heights is inferred. Extensive 
regions have been found where the surface pressure ex- 
ceeds the triple point pressure of water. 
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INVESTIGATION OF THE MARTIAN ENVIRONMENT BY INFRARED 
SPECTROSCOPY ON MARINER 9 


INTRODUCTION 

A measurement of the thermal emission spectrum yields much informa- 
tion on the environment and climate of a planet. In certain spectral invertals 
the remotely sensed radiance originates at the surface, in others the emission 
comes from atmospheric gases and aerosols. This property of the thermal 
emission permits retrieval of surface and atmospheric temperatures, surface 
and atmospheric compositions, and total amount of CO 2 in an atmospheric 
column. Other parameters, such as surface pressure and topography, and 
the wind field may then be derived. 

However, to extract this information the spectra must have adequate 
spectral resolution ( < 5 cm - l) to at least resolve details in the contours of 
gaseous absorption bands and must be recorded with adequate precision. In- 
strumentation for infrared emission spectroscopy on the Nimbus earth satel- 
lites successfully demonstrated the feasibility of this technique for meteoro- 
logical and geophysical investigations. 

Infrared emission spectroscopy on Mars has been applied from the ground 
(e.g. , Sinton and Strong, 1960) and from spacecraft (Mariner 6 and 7, Horn 
et al. 1972). However, only the orbiter Mariner 9 could explore the thermal 
emission of Mars systematically. This was the scientific objective of the in- 
frared interferometer (IRIS) experiment. The more specific goals have been 
detailed previously (Hanel et al. 1970a) and preliminary results have already 
appeared (Hanel et al. 1972a). The present paper is a further report on re- 
sults of continuing data analysis and conclusions through May 1972. 

INSTRUMENT AND CALIBRATION 

The Michelson interferometer on Mariner 9 records the spectrum be- 
tween 200 and 2000 cm“l (5 and 50 pm) with a nominal spectral resolution of 
2.4 cm‘1 (apodized) and spatial resolution of about 125 km at periapsis with 
vertical viewing. A noise-equivalent-radiance (NER) of 0. 5 x 10-7 watt cm'^ 
sr - Vcm“l has been achieved. 

The instrument is in many respects similar to interferometers flown in 
each orbit on Nimbus 3 and 4 (Hanel et al. , 1970b and 1971). However, the 
Mariner interferometer covers a wider spectral range and has a higher spectral 
resolution as well as a lower NER than the earlier Nimbus instruments. The 
Mariner interferometer has a single infrared detector (a thermistor bolometer) 
and a cesium iodide substrate for the beamsplitter coatings. The Michelson 
mirror moves with a constant velocity during the recording of an interferogram. 
Wavenumber calibration is provided by referencing all infrared spectra to the 
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0. 6929 jim line of a neon discharge tube. Other details of the design and per- 
formance of the Mariner interferometer are discussed elsewhere (Hanel et 
al. 1972b). 

Calibration spectra are periodically recorded while observing either deep 
space or an on-board blackbody. Scaling of the raw Martian spectra to the 
calibration spectra permits a display of the Martian spectra in radiometric 
units. The calibration method of the Mariner interferometer is similar to that 
of the Nimbus instruments previously described (Hanel et al. 1972c). A final 
calibration must include second order effects. Two such effects for which 
allowance was made in the calibration of the Nimbus 4 data are not present in 
the Mariner instrument. The absence in the Mariner interferometer of an 
asymmetry between the target and space ports was verified by measurements 
just before Mars encounter. The other effect in the Nimbus instrument arose 
from small temperature oscillations due to periodic loss of thermal control. 

In the Mariner IRIS the thermostat has worked well within the design range, 
controlling the instrument temperature to within a fraction of a degree. 

The only second order effect incorporated so far in the Mariner 9 calibra- 
tion procedure is necessary to account for the imperfect blackness of the warm 
calibration source, which is an aluminum plate with 30° V-shaped grooves 
painted with 3M401-C10 black velvet paint. While this paint is relatively black 
over most of the instrument spectral range, small glass beads contained in it 
give rise to emittance variations of a few percent at 1100 cm"*- and 480 cm - -*- 
which are characteristic wave-numbers of SiC>2. A correction factor was de- 
rived from laboratory reflectance measurements on a duplicate blackbody, 
from similar measurements on the same type of paint, kindly made available 
by James Aronson (private communication), and finally by comparisons of the 
warm and cold calibrations on the interferometer while in orbit around Mars. 

All three methods were in agreement and, consequently, the emissivity cor- 
rection of the warm calibration source has been applied to all spectra shown in 
this paper. 

The calibrated emission spectra are the prime product of the infrared 
spectroscopy experiment. After final calibration the Mars spectra will be 
available to the scientific community through the National Space Science Data 
Center. 

GENERAL CHARACTERISTICS OF THE SPECTRA 

Since the arrival of Mariner at Mars in November 1971 the interfer- 
ometer has produced over 20,000 spectra showing variation with Martian 
coordinates, local time, season and secular events such as the dust storm. 
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The substantial differences between spectra from the north polar hood, the 
south polar region, and mid-latitudes on Mars and from the Sahara desert 
are shown in Figure 1. 

Throughout the range from 250 to 500 cm ^ numerous rotational lines of 
water vapor are evident (Figures lb and c); the strongest lines (at 254, 278, 

303, 324, and 328 cm' 1 ) are used for the determination of atmospheric water 
vapor content. Above 1300 cm -1 weaker lines of water vapor appear also. 
Within the range from 500 to 800 cm -1 the main spectral features are due to 
numerous bands of carbon dioxide. The theoretical carbon dioxide spectrum is 
not in quantitative agreement with the observed spectra at 579, 608, 710, 730 
and 770 cm~l. Readjustment of the theoretical line parameters is expected to 
eliminate most of these discrepancies. Weak carbon dioxide bands are seen at 
828, 865, 961, and 1064 cm -1 (Figure lc). Comparison of theoretical and 
observed spectra also indicates that the isotopic ratios of C 1 -^ to C 1-8 andO 
toO 18 on Mars are close to their terrestrial values. The broad diffuse ab- 
sorption feature, Figure lc, centered at approximately 1100 cm -1 and a 
smaller feature near 480 crn - -^ are associated with dust suspended in the 
atmosphere. Another weak feature at approximately 1350 cm -1 is also tenta- 
tively identified with the suspended dust. 

The spectrum (Figure la) of the north polar area (revolution 102) closely 
resembles that of a 145K blackbody except in the C0 2 band between 600 and 700 
cm -1 . There the higher intensity indicates a warmer atmospheric layer above 
a lower boundary having approximately the condensation temperature of C0 2 . 
The south polar spectrum (Figure lb, revolution 30) has already been dis- 
cussed (Hanel et al. 1972a); the shape of the C0 2 band indicates the presence 
of a strong temperature inversion. The Martian mid-latitude spectrum (Figure 
lc, revolution 92) is compared to a Sahara spectrum (Figure le) recorded by 
the Nimbus 4 meteorological satellite. Measured over a similar spectral 
range and with nearly the same spectral and spatial resolution, these spectra 
provide a comparison between Mars and Earth. Similarities in the C0 2 band 
(Figure lc and e) indicate that temperatures decrease with altitude in both 
atmospheres except at high altitudes where the reversal of the Q-branch (667 
cm" 1 ) in the Sahara spectrum is caused by a temperature inversion due to 
heating by ozone; indication of such an inversion is lacking in the Martian 
spectrum. These results are consistent with the presence of the 1042 cm' 1 
Og band in the earth spectrum, and its absence in the Martian spectrum. Water 
vapor lines appear only weakly on Mars although strongly evident on earth. 
Shown in Figure Id is the spectral emissivity of quartz (pure Si0 2 , fractured 
solid sample). This feature may be recognized in the Sahara spectrum be- 
tween 1080 and 1250 cm -1 . The corresponding silicate dust features on Mars 
are somewhat different because of the effects of particle size, mineral com- 
position and suspension in the atmosphere, as will be discussed later. 


3 



ATMOSPHERIC WATER VAPOR 


Water vapor in the Martian atmosphere has been conclusively identified 
from rotational lines which are evident in the observed spectra. Figure lb and c. 
Preliminary analyses of the spatial and temporal changes in the water vapor 
have been made for several regions on Mars: the south polar area, mid- 
latitude regions near the sub-solar point, and over the north polar hood. South 
polar passes were analyzed from revolutions 30 (L s = 297°) and 116 (L s = 321° ) 
where L g is the planetocentric longitude of Mars, (The American Ephemeric 
and Nautical Almanac, 1971, 1972). Sub-solar passes were selected early in 
the mission on revolution 20 (L g = 293°) where the atmosphere was nearly 
isothermal, on revolution 92 (L g = 314°) during the transition period from a 
dusty to a clear atmosphere, and on revolution 174 (L g = 336°) after the at- 
mosphere had cleared substantially and the temperature lapse rate had in- 
creased to between 3 and 4 K/km. Local times were restricted to when the 
surface temperature attains its maximum value (1200-1400 hours), thus maxi- 
mizing the thermal contrast between surface and atmosphere; this in turn en- 
hances the H 2 0 lines in the observed spectra. For each water vapor determina- 
tion ten to fifteen spectra are averaged to improve the signal-to-noise ratio. 

The total H 2 0 content is determined from a quantitative comparison of 
observed and synthesized spectra. The synthetic radiances were calculated by 
a numerical integration of a solution of the radiative transfer equation. Tem- 
peratures obtained from the inversion of the 667 cm~l CO 2 band, as discussed 
below, and a constant H 2 0 mixing ratio were used. The latter is assumed based 
on the evidence of strong vertical and horizontal mixing indicated by the ex- 
istence of the dust storm and by the wind fields discussed below. Synthesized 
and observed spectra for the south polar area are compared in Figure 2; the 
former are offset vertically for clarity. Synthesized spectra with 10 to 20 
precipitable micrometers of water vapor give best agreement with the ob- 
servations. Similar comparisons for the sub-solar point cases also yield 
abundances in the range of 10 to 20 precipitable micrometers. Over the north 
polar hood water vapor was not identified, presumably due to the inability of 
the atmosphere to hold significant amounts of water vapor at the low tempera- 
tures. It should be emphasized that estimates of the water vapor content are 
extremely sensitive to the surface pressure, the temperature difference be- 
tween the surface and the lower atmospheric layers, and to the temperature 
gradient in the lower portion of the atmosphere. With refinements in these 
parameters the estimated amount of water vapor is also subject to adjustment. 

Water vapor results from both IRIS spectra and ground-based observa- 
tions are summarized in Figure 3. The upper portion of the figure presents 
data from sub-solar point latitudes (-20 to -30°) and the lower portion from 
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the south polar region. The values determined from IRIS are indicated by 
vertical bars. The cross-hatched areas represent the range of total H 2 0 con- 
tent determined from ground-based observations during the 1969 opposition by 
Barker et al. , (1970). The horizontal dashed lines represent mean values 
from the ground-based observations made by Tull (1972) during the same time 
period of the IRIS measurements. All ground-based observations represent 
averages over the entire disk of Mars. In contrast to the above measurements, 
the results from the 1.38 pm water vapor band experiment on Mars 3 (Moroz 
and Ksanfomaliti, 1972) indicate substantially lower amounts of water vapor. 

From these results it may be concluded that the average abundance of 
atmospheric water vapor decreased slightly as the entrained dust settled but 
then increased after atmospheric clearing during the late southern summer, 

L g = 350 to 360°. The latitudinal and seasonal changes of atmospheric water 
vapor on Mars are often correlated with changes in the polar caps (Barker 
et al. , 1970; Schorn, 1971). However, in 1971 the water vapor content did not 
decrease towards the equator from the receding south polar cap as would be 
expected if it were the source of vapor. Latitudinal gradients were noted from 
the north polar cap during the 1969 opposition (Tull, 1970). Perhaps most 
significant, however, is the fact that the average abundance of water vapor for 
the entire planet is lower than in previous oppositions. This lower abundance 
could be due to an unusually large amount of water trapped in the north polar 
cap; some evidence of anomalous behavior during the northern autumn is evi- 
dent from the extension of the north polar hood to lower latitudes than in pre- 
vious corresponding seasons. 

Alternatively, however, the unusually low water vapor content may be 
associated with adsorption on dust particles. Experimental and theoretical 
investigations by Davis (1969), Pollack et al. (1970), and Fanale and Cannon 
(1971) have indicated that adsorption of water vapor by the Martian regolith 
may be a significant mechanism in governing the vapor abundance. In general, 
gas is adsorbed at low temperatures and desorbed at higher temperatures. It 
is suggested that the suspended dust provided additional surface areas to per- 
mit more efficient adsorption and that much of the water then settled with the 
dust. With increasing temperatures at the surface some of the water vapor 
may return to the atmosphere. 

As more and more refined water vapor analyses are available a more 
quantitative investigation of the relative importance of adsorption and of the 
polar cold trap in these processes will be made. It should also be noted that 
water vapor amounts of 10 to 20 precipitable micrometers, in conjunction with 
the inferred temperature profiles, allow condensation in the upper atmosphere. 
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ATMOSPHERIC DUST COMPOSITION 


South-polar and mid-latitude spectra recorded during the first 100 orbits 
show strong features whose positions and relative strengths are characteristic 
of silicate minerals. The broad features near 480 cm -1 and between 900 and 
1200 cm -1 - which appear in emission in the south-polar spectra, and in absorp- 
tion over the low latitudes where the surface is warmer than the atmosphere, 
indicate the presence of particles suspended in the atmosphere. It has been 
pointed out that the mineralogical characteristics of silicates (primarily the 
silicon-oxygen grouping) determine the spectral position of the absorption and 
transmission maxima, (Lyon, 1964; Hovis and Callahan, 1966; Aronson et al. 
1967; Conel, 1969; Salisbury, 1970). From such data it may be possible to 
infer the dust composition. However, the matter is complicated by effects of 
particle size and the possibility that more than one mineral group may be con- 
tributing to the recorded spectrum. 

The wavenumber of maximum absorption in the interval between 900 and 
1250 cm -1 may be obtained in several ways. One is to determine the minimum 
in the brightness temperature of the mid-latitude spectra and the corresponding 
maximum in the south polar spectra; samples of both are shown in Figure 4. 
The brightness temperature is better suited for such a determination than the 
intensity because the nonlinear effect of the Planck function has been removed. 
A value between 1070 and 1100 cm - l may be deduced from the spectra. 

A second method, based on the standard deviation of a set of spectra, 
yields a more precise value. While the shallow minimum of the brightness 
spectrum could conceivably be affected by unrecognized instrumental effects, 
the standard deviation depends much less on the absolute calibration. Shown 
in Figure 5 is the same mid-latitude average spectrum (Rev. 92) displayed in 
Figure 4a, together with the standard deviation cr (T). The ensemble of spec- 
tra was chosen between 1200 and 1400 hours local Mars time. Over this 
period, the surface temperature is near its diurnal maximum and is expected 
to change very little between individual samples. The lower atmosphere, 
however, is still warming as indicated by the maxima in the standard deviation 
near 600 and 730 cm“l. At these wavenumbers atmospheric emission origi- 
nates primarily between the 3 mb level and the surface. A low standard devia- 
tion in the center of the 667 cm -1 - CO2 band indicates stable temperatures at 
higher levels. The strong maximum at 1090 cm' 1 is caused in part by tem- 
perature fluctuations in the lower atmosphere but also by the local variability 
of the dust among spectra in the ensemble. In either case the maximum varia- 
tion is expected to occur at the most absorbing wavenumber. This occurs at 
1090 cm -1 -, in agreement with the results obtained from the brightness 
temperatures. 
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It has been pointed out by Conel (1969) that the transmission maximum 
(Christiansen frequency) is also a good indicator of bulk chemical composition. 
For silicates, the Christiansen frequency shifts from 1360 crn'l for quartz 
towards lower wavenumbers with the inclusion of other elements into the 
lattice structure. From the maximum of the brightness temperature in the 
mid-latitude spectrum and from the minimum in the standard deviation (Figure 
5) a Christiansen frequency between 1280 and 1300 cm -1 is inferred. Again, 
both methods are consistent. 

Although the spectral position of the absorption maximum and minimum of 
the Martian dust cloud seems to be well defined, the interpretation in terms of 
the silica content is still difficult. The difficulties arise in part from the com- 
plexity of theoretical treatments of a cloud of scattering and absorbing particles 
of non-spherical shape and in part from the lack of knowledge of the complex 
refractive index of most minerals in the infrared. Although some effort to- 
wards a better theoretical understanding have been made by Aronson et al. 

(1967) and Conel (1969) the interpretation of the IRIS spectra in terms of the 
Si0 2 content must still be based on an empirical comparison with laboratory 
transmission spectra of mineral dust. 

One extensive set of transmission measurements of mineral dust has been 
prepared by Lyon (1964) who imbedded mineral powders in KBr pellets. With 
such preparation only the absorption maxima are reliable, since the minima 
representing the Christiansen frequency are shifted considerably by the dif- 
ferences in refractive index between air or CO£ (~ 1) and KBr (~ 1.5). Com- 
parison of the peak absorption of the Martian dust as measured by IRIS 
(~ 1085 cm -1 ) with the transmission spectra of Lyon indicates a reasonable 
agreement with dacite (68.7%SiC>2> USNM 82, Lassen Peak, Cal.) as well as 
with quartz basalt (57. 25% Si0 2 , USNM 101, Snag Lake Cinder Cone, Lassen 
Co. , Cal. ) but poor agreement with more basic materials. 

Comparison of the Christiansen frequency measured by IRIS (~1290 cm -1 ) 
with Conel' s laboratory data (Conel 1969 Figure 5) shows good agreement with 
the alkali-feldspar group (Si0 2 ~ 65%) as well as with andesine (Si0 2 ~57%) 
and oligoclase (SiC>2 ~ 62%) of the plagioclase group, but no agreement with 
pure SiC>2, or with micas or other minerals of lower Si0 2 content. Comparison 
of the IRIS-determined Christiansen frequency with a diagram kindly made 
available by J. Salisbury (private communication) indicates also a Si0 2 content 
by weight of approximately 60%. 

From these comparisons the Si0 2 content of the Martian dust may be 
estimated as 60 ±10%. 
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GEOCHEMICAL IMPLICATION 


Before drawing any geochemical conclusions from the chemical com- 
position of the Martian aerosol, it is important to determine if the dust, either 
suspended or on the surface, is chemically representative of the Martian sur- 
face rocks, and specifically if the dust could be higher in Si 02 than the rocks 
from which it was derived. 

A possible parallel to the Martian situation may be provided by ter- 
restrial wind-deposited sediments (loess), found in large areas of North 

America, Europe, and Asia. Pettijohn (1957, p. 378) lists chemical analyses 
of 5 loess samples from the United States, China, and Belgium, with the fol- 
lowing weight percentages of Si0 2 : 59.3, 61.0, 64.6, 72.8, and 74.5. These 
values are somewhat more siliceous than the average continential crust, which 
is about 60% Si0 2 according to studies by Ronov and Yaroshevsky (1969). How- 
ever, it is apparent that no gross chemical fractionation occurred in these 
sediments. 


A second aspect of the question involves possible chemical fractionation 
by grain size of wind- deposited sediments. As Pettijohn (1957, p. 101) shows, 
chemical composition of terrestrial sediments (including water-deposited ones) 
is strongly dependent on grain size. Fine-grained sediments tend to be richer 
m clay minerals (chiefly kaolinite, about 46% Si0 2 ), and coarse-grained ones 
richer in quartz (pure Si0 2 ) ; thus fine-grained sediments are generally less 
siliceous than coarser-grained ones. For' example, terrestrial loess has 
grain sizes largely between 30 and 60 micrometers (Pettijohn, 1957, p. 378), 
compared' with the average grain size of less than 5 micrometers for clay 
(Krumbein and Sloss, 1951, p. 71), a major constituent of terrestrial sedi- 
ments as a group. Evidence by IRIS and other experiments (Moroz and 
Ksanfomaliti, 1972) indicates that the dust particle size is a few micrometers 
It is therefore expected that the dust is probably somewhat less siliceous than" 
the rocks from which it is derived. 

Fr ° m the above evid ence, it appears that the present estimate of 60 ±10% 
i 2 f° r the Martian dust composition is probably fairly representative of the 
mean composition of the Martian crust. This being so, it supports the con- 
c us ion that geochemical differentiation of Mars has occurred (Hanel et al. , 
1972a). Comparison of various differentiated and undifferentiated bodies is 
shown m Table 1. The conclusion that Mars is differentiated is further sup- 
ported by the existence of volcanoes and lava fields discovered by the television 
experiment on Mariner 9 (McCauley et al. , 1972). Many of the lava flows, 
such as those in the Nix Olympica-Tharsis region, appear similar in 
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morphology and scale to basaltic lava flows as seen on terrestrial orbital photo- 
graphs (Lowman and Tiedemann, 1971). 


Table 1 


Comparative Si0 2 Content of Terrestrial, Lunar, 
and Martian Crusts and Chondrites 


Earth 


Moon 


Mars Chondrites 


Continental Oceanic Maria 2 Highlands 
SiO? (wt.%) 60.2 48.7 39.9 47.1 60±10 


46.8 


1 Wyllie, 1971. 

2 Lowman, 1972; average of 22 crystalline rocks from Apollo 11 and 12 sites. 

3 Average of 8 fragmental rocks from Apollo 14 site (Lowman, 1972) and 2 soil samples from Apollo 15 
Apennine Front site (Apollo 15 Preliminary Science Report, 1972). 

4 Richardton chrondrite (analyst, H. Wiik), cited as typical by Wood (1968); Si0 2 given relative to stony 
material; makes up 34.3% (wt.) of total meteorite including metallic and sulfide phases. 


SURFACE PRESSURE 

The atmospheric pressure at the Martian surface, and from this the 
topography of the planet, is determined by evaluating the opacity of the atmos- 
phere in the wing of the 667 cm’'- band of CO^. 

The upwelling radiance I„ observed at frequency v is given by 

/ 0 T 

By [T (In p)] ■■ " ~ d In p (1 

v 31np 

where e„ is the emissivity of the surface, B„ (T g ) is the Planck function for 
radiance emitted by the surface at temperature T g , t uS is the transmittance 
of the atmosphere from the surface to the spacecraft, and the integral is the 
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contribution of the atmosphere. By rewriting the integral in terms of an 
average over transmittance \ , (1) may be rewritten as 

l v - ( B „) 

T "S - € v B p (T s > - (B„) < 2 > 

This expression is used iteratively to evaluate t„ s . l y is the directly measured 
radiance, T s is evaluated from intensities in the window region (1300 cm - -*-) of 
the infrared spectrum, and is evaluated just outside of the 667 cm -1 band 
and is assumed to remain constant within the edge of the band. It is assumed 
that the atmospheric thermal structure does not change significantly between 
adjacent IRIS measurements, so that the temperature profile from a given 
measurement may be used in estimating T „ s for the following measurement. 

The surface pressure is obtained from t vs by assuming a C0 2 mixing ratio of 
unity and employing a transmittance model based on direct line-by-line integra- 
tion techniques. & 


Topography of the Tharsis area in 1 km increments is shown in Figure 6. 
Altitudes are referred to a 6. 1 mb pressure surface, and are evaluated from 
the surface pressure using a 10 km scale height. The contours are alternately 
dotted and solid for clarity. The IRIS instrument did not view North Spot (11° , 
104°), so it is not shown; the instrument viewed only the slope of Middle Spot 
(0 , 113 ), so its height could not be determined; a complete spectrum was 
obtained in the caldera of South Spot (-9°, 120°), from which the altitude of the 
summit crater is found to be 14 km. Considering the 8 km contour as the base 
of the cone, the average slope of the structure is approximately 3° . This low 
slope suggests that the crater was formed by very fluid lavas, which is con- 
sistent with the interpretation of these features as shield volcanoes (McCauley 
et al. , 1972). A further feature of interest is the valley which proceeds east- 
ward from (-8°, 98°), which corresponds to the beginning of the great equa- 
torial rift canyon (McCauley et al. , 1972). 

Previous spacecraft measurements have shown that Hellas is a sub- 
stantial depressed area (Herr et al. , 1970; Barth and Hord, 1971). Figure 7 
shows the topography of Hellas basin as obtained from IRIS data. Considering 
the feature to be defined by the 0 km contour, it is found to be roughly 1800 km 
in its north-south dimension, and 2200 km in its east-west dimension. The 
eastern slope is noticeably less steep than that on the west. Comparison of 
these contours with the TV pictures of the basin shows that the transition from 
rugged terrain to the relatively featureless terrain deep within Hellas occurs 
at about the -1km (6. 7 mb) level. 
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Areas in which the surface pressure exceeds 6. 1 mb are of biological in- 
terest, for under temperature conditions which frequently occur on Mars, 
liquid water can exist in equilibrium with its vapor. Further preliminary 
topographic mapping between -65° and +20° latitude discloses that the surface 
pressure during southern summer (320°< L s < 350°) exceeds 6.1 mb in fairly 
large areas of Argyre, western Margaritifer Sinus and Isidis Regio, in addi- 
tion to Hellas. 

Preliminary comparison of IRIS topography with 2 pm CO 2 (Herr et al. 
1970, Moroz and Ksanfomaliti, 1972) and ultraviolet results (Barth and Hord, 
1971) generally show good agreement for altitudes above approximately 0 km. 
•However, results of the above experiments are systematically higher than 
IRIS results at lower altitudes, for example, in Hellas. Possible causes of the 
discrepancies, e.g. , the presence of residual fine dust in the atmosphere, are 
being investigated. 

Comparison of IRIS topographic altitudes with occultation (Kliore et al. 
this issue) and with radar (Downs et al. 1971; Pettengill et al. 1971) results 
shows the following. Below the 4 km level, IRIS topographic data lie sys- 
tematically 1 km below occultation results. At increasingly high altitudes, 
however, IRIS data fall progressively lower than occultation data, the dif- 
ference reaching 3 km at an occultation altitude of 9 km. A similar trend is 
observed in comparison with radar data, in addition to an overall offset of 
about +3 km at the radar zero of altitude. Aside from these systematic dif- 
ferences, the scatter in the IRIS data relative to either occultation or to radar 
measurements is about ±1 km. At present, the CO 2 transmittance functions 
and irregularities in the Martian gravitational potential field (Lorell et al. , 
1972a, 1972b) are being investigated as possible sources of the systematic 
differences. Further possibly significant sources of error are pressure varia- 
tions due to cyclonic activity, to diurnal tidal effects, and particularly for the 
extended mission, to seasonal effects associated with polar cap changes 
(Leighton and Murray, 1966). 

ATMOSPHERIC AND SURFACE TEMPERATURE 

Another objective of the IRIS experiment is to derive information on the 
thermal structure of the Martian atmosphere, using the 667 cm-1 C0 2 absorp- 
tion bands. Data from the strongly absorbing band center provide information 
on the temperature of the upper atmosphere while data from the less strongly 
absorbing portions of the bands provide information on the middle and lower 
atmosphere. A knowledge of the spectral transmittance of the atmosphere 
above each atmospheric level is required. Transmittance models based on 
direct line-by-line integration techniques have been employed for this purpose. 
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In the studies carried out thus far, the Martian atmosphere has been assumed 
to be entirely CO 2 - The techniques employed for inverting the integral equa- 
tion of radiative transfer (1) to obtain T(p) are similar to those developed for 
remote sounding of the terrestrial atmosphere (see Conrath and Revah, 1972 
for a review). 

Preliminary estimates of temperature profiles have been made for almost 
all usable spectra thus far obtained. It is anticipated that future estimates will 
incorporate some revisions in the atmospheric transmittance models, im- 
provements in the calibration of the spectra, and possibly some modifications 
in the computational techniques. Therefore, the results presented here should 
still be regarded as preliminary. However, it is felt that the gross properties 
of the inferred atmospheric temperature field will not be modified by sub- 
sequent improvements in the analyses. 

Examples of three individual temperature profiles are shown in Figure 8. 
All three profiles were acquired at approximately the same latitude (-30°) and 
local time (1400 hours). The profile from revolution 20 is typical of those ob- 
tained in the early part of the mission when the atmosphere was heavily dust- 
laden. The lower part of the profile is regarded as uncertain since possible 
additional opacity in the 667 cm"-'- band due to dust was not included in the 
analysis. As would be expected, profiles obtained during this period do not 
agree well with theoretical predictions based on dust-free atmosphere (Gierasch 
and Goody, 1968; Prabhakara and Hogan, 1965; Ohring and Mariano, 1965). 
However, reasonable agreement is obtained with recent calculations of Gierasch 
and Goody (1972) which include grey atmospheric absorption of solar energy. 

As the dust began to dissipate, the apparent lapse rate in the lower atmosphere 
increased while the upper atmosphere cooled as indicated by the example from 
revolution 92. By revolution 174, further cooling had occurred at all levels. 
However, in all cases the lapse rate remained sub- adiabatic. 

Initial analyses of the temperature field have focused on obtaining the 
mean diurnal behavior and its change as the dust dissipated. Data from se- 
lected groups of orbits were sorted into sets spanning 10 degrees in latitude 
and one hour in Martian local time. The average was taken for each of these 
sets and used to construct isotherms in a latitude-local time plane. In this 
way, the mean dirunal behavior of the temperature field could be estimated for 
various atmospheric levels. 

Examples from three different groups of orbits are shown in Figures 9, 
10, and 11. In each case only two atmospheric levels are shown, the 0.3 mb 
level (~30 km) and the 2 mb level (~10 km). While large gaps in the data 
exist, the gross behavior of the temperature field is apparent. The data shown 
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in Figure 9 were obtained between orbits 1 and 85 when a large amount of dust 
was present in the atmosphere. Temperature maxima are seen at both levels 
near latitude -60° late in the afternoon. Between orbits 85 and 120 (Figure 10) 
atmospheric clearing was evident, and the temperature maxima decreased in 
magnitude and shifted toward the equator and earlier local times. By orbits 
161 - 186 (Figure 11) further cooling had occurred and the maximum at the 2 
mb level had shifted to approximately the latitude of the sub-solar point (-7°). 
The observed diurnal variations of 15 to 30K are substantially larger than 
theoretical predictions under both dusty and clear conditions (Gierasch and 
Goody, 1972). 

A similar analysis of the diurnal behavior of the brightness temperatures 
at 1300 cm _ l was carried out. Near this wavenumber the mineral dust seems 
to be most transparent (Christiansen frequency) and therefore the measured 
brightness temperature is the best estimate of the temperature of the surface. 
As shown in Figure 12, the surface temperatures are found to behave quite 
differently from the atmospheric temperatures. The maximum occurred near 
the latitude of the sub-solar point at about 1300 hours local time both during 
the period the atmosphere was dust filled and after the clearing was essentially 
complete. 

WIND FIELDS 

Preliminary estimates of wind fields consistent with the observed mean 
atmospheric temperature fields have been made. Analytical models of the 
longitudinally averaged temperature fields were constructed for use in fluid 
dynamic equations from measured atmospheric temperatures obtained in 
orbits 1 to 85 and 161 to 186. The model temperatures consist of an axially 
symmetric component and a sinusoidally varying diurnal component. The lapse 
rate and the location of temperature maxima in altitude, latitude and local time 
have been incorporated into the models so that the gross temperature gradients, 
as determined from the spectra, have been retained. The temperatures have 
been extrapolated into the regions of local time and latitude where data are 
missing. 

The approximate hydrodynamic equations used to describe atmospheric 
motion in terms of latitude 0, longitude X , and pressure p, are (see e.g. , 
Haltiner, 1971). 


dV 

dt 


= -V$ 


- fk x V + F u 
rl 


(3) 
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The hydrostatic approximation has been assumed, and the vertical component 
of the Coriolis acceleration has been neglected. V is the horizontal velocity, 

V is the two dimensional gradient operator on a constant pressure surface, k 
is a unit vector in the direction of the local vertical, and f = 2S2 sin ip is the 
Coriolis parameter for a planetary rotational angular velocity Q . Fh is the 
horizontal frictional force, $ is the geopotential function, and in the coordinate 
system employed, the parameter co = dp/dt replaces the vertical velocity. 


The geopotential function can be obtained from (5) in the form 

P 

$(«/>,*., P, t) = -R J T (ip, p, t) din p + $ g (ip, X) (6) 

P s (ip,h » t) 

where <j> g is the geopotential at the planetary surface. In general, $ s and p s 
will vary with-latitude and longitude because of changes in topography, and p s 
will also change due to atmospheric motions. Because of the large elevation 
differences found on Mars, the effects of topography on the wind field may be 
considerable (Gierasch and Sagan, 1971; Sagan et al. , 1971; Blumsack, 1971). 
In addition, the large observed diurnal variations in the atmospheric tempera- 
tures can result in significant tidal variations in p s . However, in this pre- 
liminary study $ s and p s are assumed to be independent of latitude and longi- 
tude. When the mean temperature field is then employed in (6), a geopotential 
function is obtained which is dependent only on latitude, pressure, and local 
time. Further approximations include neglect of the nonlinear transport terms 
in (3) and the assumption that F h = -T “ly. The damping time r was chosen 
as 4 days, corresponding approximately to an eddy transport coefficient of 2 x 
107 cm^/ sec. With these assumptions, the greatly simplified momentum 
equation is 

P 

R f — dp - f k x V - T -1 V (7) 

P 

P s 

Numerical integration of (7) was carried out over times long compared to 
t , resulting in a steady state solution independent of the arbitrary initial 
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conditions employed. The wind field at one pressure scale height resulting 
from the temperature model based on the mean temperatures from orbits 1-85 
is shown in Figure 13. Winds obtained in the equatorial region are highly un- 
certain due to the vanishing of the Coriolis terms and dominance of the damping 
terms. The most significant feature of the wind field is the strong diurnal tidal 
component. It is of interest to compare this wind field with that obtained from 
(7) with 3 V/ 9 t set equal to zero. The resulting equation is the thermal wind 
equation with an additional damping term which is significant only near the 
equator. The thermal wind field calculated using this equation is shown in 
Figure 14. Comparison of the two wind fields indicates that the thermal winds 
are a rather poor approximation, and there is considerable cross-isobaric flow 
due to the tidal nature of the wind regime. 

Wind fields were also calculated from (7) using the mean temperature 
field from orbits 161-186. A strong tidal component was again obtained. The 
stronger meridonial temperature gradients resulted in somewhat larger wind 
velocities in this case. 

The highly simplified dynamical model considered here can be expected 
to give only the gross behavior of the mean wind field. More sophisticated 
models which include the frictional interaction of the atmosphere and surface, 
non-linear transport terms, tidal variations in surface pressure, and topography 
are currently being developed. 

BIOLOGICAL IMPLICATIONS 

Mariner 9 has produced data suggesting that the prospects for life on 
Mars, deemed negligible since Mariner 4, are significantly increased. IRIS 
has shown that differentiation has occurred on Mars. Dynamic events such as 
dust storms, high winds, hazes or clouds, and polar cap transformations have 
also been recorded. Mariner 9 television data (McCauley et al. , 1972) have 
indicated that faulting, volcanism, wind and water erosion have been active 
forces on the planet. Many of these processes are conducive to the evolution 
of life. For example, differentiation of the planet makes lighter elements 
more readily available for evolutionary processes at the surface. Dynamic 
processes ensure the exposure and mixing of possible life precursors and 
nutrients. Further, they create a diversity of habitats essential for divergent 
evolutionary processes. 

Trace atmospheric constituents, which might possibly imply surface 
metabolic processes in accordance with the theory of Hitchcock and Lovelock 
(1967), and Dayhoff et al. (1967), have so far not been identified in IRIS spectra. 
However, the observed atmospheric composition, atmospheric pressure, 


15 



temperature, and surface conditions do not preclude the existence of an in- 
dependent life on Mars. The general lack of significant amounts of water has 
been considered the major life limiting factor (Mariner-Mars 1969: A Pre- 
liminary Report, 1969). Since liquid water may currently be lacking on the 
planet surface and only traces of water vapor exist in the atmosphere, Martian 
life would have had to adapt to extreme aridity or have become restricted to 
limited areas where liquid water may be available at least transitorily. Topo- 
graphic results discussed above indicate that liquid water could exist in exten- 
sive areas on the surface of Mars if adequate sources of juvenile water were 
present; such areas then represent the most likely places to explore the possible 
existence of life. 

It is generally agreed that, if life has evolved on Mars, both phototrophic 
and heterotrophic organisms are required to harness the sun's energy and to 
maintain nutrient cycles (Morowitz, 1968). Assuming the existence of microbial 
life, the dynamic surface forces indicated by Mariner 9 suggest that hetero- 
trophic microbial forms (and chemoautotrophs) would be widely distributed. 
Microbial phototrophs, on the other hand, may have been selected against and, 
if extant, would probably be restricted to microhabitats where they are pro- 
tected from the shifting dust which would otherwise exclude the sun's energy. 

Observed Martian conditions do not preclude the survival of some known 
earth life forms although their growth under these conditions is speculative. 
Many organisms, in fact, are known to survive conditions similar to those on 
Mars with respect to temperature, water content, and pressure. Further, 
terrestrial microbes can readily be cited which grow in saline water, acidic 
and alkaline solutions, at freezing temperatures, or anaerobically under carbon 
dioxide. These considerations, coupled with the dynamic mechanisms available 
on the planet for spreading contaminants, emphasize the necessity for steriliza- 
tion of Mars landers. 

SUMMARY 

Preliminary analysis of Mariner 9 IRIS data through May 1972 has thus 
far revealed spectral features due to CC>2» water vapor, and silicate dust sus- 
pended in the atmosphere. A number of significant results have been obtained 
and conclusions have been reached. 

Water vapor abundances of 10-20 precipitable micrometers were found 
both during the dust storm and as the dust dissipated, in general agreement 
with ground-based measurements obtained during the same period. However, 
ground-based measurements obtained in previous years during similar Martian 
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seasonal conditions have indicated larger amounts of water vapor. No sig- 
nificant latitudinal variation in water vapor was observed from the south pole 
to the equator. 

A study of spectral features due to the dust indicates a Si02 content of 
60% ±10%. Arguments have been given which indicate the surface material 
should be at least as acidic as the atmospheric dust. These results support 
the conclusion that substantial geochemical differentiation has occurred on 
Mars. 


Surface pressures have been estimated using measurements in the wing 
of the 667 cm -1 C0 2 band. Preliminary results demonstrate the usefulness of 
this technique for obtaining surface topography over large portions of the planet. 
Topographic variations of the order of 2 pressure scale heights have been 
found, and large regions exist for which the pressure lies above the triple point 
pressure for water. 

Atmospheric temperatures, also derived from the 667 cm ^ COg bands, 
always showed subadiabatic lapse rates as well as unexpectedly large diurnal 
variations of 15 to 30K. The maximum temperatures occurred late in the after- 
noon near latitude -60° during the dust storm, but moved toward the subsolar 
latitude and shortly after local noon as the dust cleared. Unlike the atmos- 
pheric temperatures, the surface temperatures displayed a maximum shortly 
after local noon and near the subsolar point at all times. 

Winds derived from the temperature fields were found to possess strong 
diurnal tidal components. Comparison of thermal winds with winds obtained 
from the linearized hydrodynamic equations indicates the existence of signifi- 
cant cross-isobaric flow. 

At present, no results indicative of biological activity have been obtained. 
However, conditions for life on Mars have been found to be less harsh than pre- 
viously thought. These findings validate the continuing efforts in the search for 
life on Mars. Further, the possibility that terrestrial microorganisms may 
survive, if not actually grow, on Mars, emphasizes the importance of sterili- 
zation of future landers. 
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Figure 2. The rotational spectrum of water vapor in the Martian 
atmosphere as seen in emission over the south polar cap 
(revolutions 30 and 116). Shown above the measured 
spectra are synthesized spectra containing 10 precipitable 
micrometers of water vapor. 
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Figure 3. Summary of the water vapor content of the Martian 
atmosphere, as determined from IRIS spectra in 
1971/1972 for the southern hemisphere summer 
season, and from ground-based observations during 
the 1969 and 1971 oppositions. The upper portion 
of the figure is for mid-latitudes, the lower portion 
is for the south polar region. 
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Figure 4. Brightness temperatures calculated from spectra recorded on revolutions 30 and 92. The mid-latitude spectrum shows a minimum due to silicate dust near 
1085 cm* 1 . The south polar spectrum shows a maximum at approximately the same wavenumber. 
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Figure 5. a) Brightness temperature average of 16 spectra recorded on revolution 92. b) Standard deviation of the 
brightness temperatures within the ensemble. The maximum in the brightness temperature near 1300 cm'1 is 
caused by a transmission maximum in the dust cloud. The maximum in the standard deviation occurs at the 
most absorbing wavenumber of the dust. 
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Figure 6* Topography of Amazonis-Tharsis region in kilometers above 6.1 mb pressure surface as determined from revolutions 142 to 160, North Spot 
(11°, 104°) was not viewed by the instrument; the flank of Middle Spot (0°, 113°) was observed; the 14 km figure at (-9°, 120^ 
represents the altitude of the caldera atop South Spot. Alternate contours dotted for clarity- 
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Figure 7. Topography of Hellas region in kilometers above 6.1 mb pressure surface as determined from revolutions 110 to 132. Alternate contours dotted 
for clarity. 
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Figure 8. Temperature profiles obtained by inversion of measurements in the 667 cm" 1 CO2 absorption band. The profiles were obtained at approximately -30° 
latitude and 1400 hours Martian local time during the revolutions indicated. During revolution 20 the atmosphere was heavily dust laden. 

By revolution 92 the dust was dissipating; clearing was well advanced by revolution 174. 
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Figure 9. Dependence on latitude and Martian local time of atmospheric temperatures at the 2 mb level 10 
km) and at the 0.3 mb level («"*■» 30 km). The temperatures represent averages over data obtained on 
revolutions 1 to 85 during the dust storm. 
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Figure 10. Same as Fig. 9, except that data used were obtained as the atmospheric dust was dissipating on 
revolutions 85 to 120. 
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Figure 12. Dependence of surface temperatures on latitude and Martian local time. Averages obtained during the 
dust storm (revolutions 1 to 85) are compared with averages obtained after atmospheric clearing had 
occurred (revolutions 161-186). 
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Figure 13. Wind fields at one scale height generated from a temperature model of the atmosphere based on data obtained during revolutions 1 
through 85. The wind velocities were calculated using the assumption of hydrostatic equilibrium and the linearized momentum 
equations, with the frictional forces represented by a linear damping term. The lower boundary was a geopotential surface on which 
a constant pressure was assumed. The wind direction is away from the dots. A vector of length equal to the horizontal or vertical 
separation of the dots represents a wind speed of 50 m/sec. 
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Figure 14. Thermal winds generated for same conditions as those of Fig. 13. The thermal winds are calculated using a balance between the 
Coriolis, pressure, and linear damping terms. 



